Final stare interaction enhancement effect on the near threshold pp system in 

S ± — > ppn^ decay 
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We discuss the low-mass enhancement effect in the baryon-antibaryon invariant mass in three- 
body baryonic B decays using final state interactions in the framework of Regge theory. We show 
that the rescattering between baryonic pair can reproduce the observed mass spectrum. 
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I. INTRODUCTION 

The mesonic B decays have been intensively studied principally to test Standard Model and to investigate the CP 
violation mechanism as predict by the Cabibbo-Kobayashi-Maskawa model. On the other hand, the large mass of the 
6-quark allows B meson to decay into a pair baryon-antibaryon too. The first measurements of baryonic B decays 
were reported by ARGUS collaboration [1], at the end of '80 years, and stimulated extensive theoretical studies. 
Interest in this area was revitalized in the last years thanks to the new measurements collected by CLEO, BELLE 
and BaBar. The phenomenology of baryonic B decays is rich and diversified (for recent reviews see Ref.[2-4]): in 
this paper we will refer to three-body B decays with a baryon-antibaryon pair and a meson in the final state. 

There is a common and unique feature for baryon-antibaryon-meson B decays, to which this paper is devoted, 
that is the observed peak, near to the threshold area, in the invariant mass spectrum of baryon-antibaryon system. 
The first experimental observation of this enhancement came from Belle collaboration studying the proton-antiproton 
system for the decays B — > ppK [5] and B — > ppD^ [6]. Very recently BaBar collaboration confirmed the threshold 
enhancement in the same channels [7, 8] with a very high statistics. Next the same peak has been found by BELLE 
studying the channels B — » ppw and B — > ppK* [9]: these results are reported in Fig. 1. 

The low-mass enhancement effect in three-body baryonic B decays indicates, thus, that there is a favorable experi- 
mental configuration, the baryon-antibaryon pair with low invariant mass accompanied by a fast recoil meson, that is 
preferred respect to the other. Therefore, the main consequence of the threshold effect is that the two-body baryonic 
decays, which have an invariant mass fixed to tub, should be disfavored respect to three-body decays. In fact it has 
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FIG. 1: Experimental results for B + — > pfm + and B + — > ppK* + respectively on the left and right panel. The shared distribution 
is from the phase-space MC simulation with area normalized to signal yield [9]. In both cases is evident the sharp peak near 
to m P p « 2m p ps 2 GeV. 
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been well established experimentally that 

B(B — ► BB'M) > B(B BB') ; (1) 

typically B{B -> BB'M) w ©(lO" 6 ) and B(B -► BB') « e>(10~ 7 ) [2-4]. In other words, thanks to the threshold 
effect, in _B — > BB'M decays the effective mass of baryon-antibaryon is reduced, respect to the baryonic two-body 
decays, as the emitted meson carries away much energy. The hierarchy on the Branching ratio in (1) is thus an 
indirect support to the enhancement effect. 

Several theoretical speculations have been proposed to interpret the anomalous observed peak at the threshold. 
The first discussion was given by Hou and Soni using a simple pole model [16], and subsequently this mechanism 
has been addressed in the QCD naive-factorization approach [17, 18]. Various interpretations include models with 
baryon-antibaryon bound state based on the old idea of Fermi and Yang [19], while exotic interpretations are given in 
terms of glueball intermediate states and fragmentation picture as in Ref. [20]. Furthermore, the problem has been 
studied also in the framework of perturbative QCD: in this scheme the enhancement is led back to the asymptotic 
behaviour of baryonic form factors [21]. 

In this paper we will follow a different philosophy. The threshold enhancement effect in proton-antiproton invariant 
mass, is not an exclusive prerogative of B meson decay but it has been observed also in other meson decays, with a very 
similar behavior: in the J/ip — > jpp decay by BES collaboration [10] and from CLEO measurements in T(IS') — * jpp 
decay [11]. Furthermore, recent observations made at the LEAR at CERN, studying the proton-antiproton production 
in the cross section of the reaction e + e~ — > pp [12], showed a near threshold structure similar to the previous one. 
Moreover the characteristic near threshold enhancement is not only present in the proton-antiproton system but also 
in the systems that contain hyperons as Ap and Ap, as observed by BELLE in B° — > pAir~ decay [13] and by BES in 
J/ip — > pAK~ + c.c. decay [14], and in the system AA as observed by BELLE in B + — > AAK + decay [15]. 

As one can see from the list of experimental data, the enhancement effect is universal and it is present with similar 
feature in different contexts: B decays as well as J/ip and T decays, it is found in the pp system as well as in Ap 
system etc. . .. Therefore these experimental results suggest that: 

the low-energy enhancement in baryon-antibaryon mass spectrum is mostly related with the dynamics of 
the baryonic pair and is weakly depending by the production decay vertex. 

For this reason we consider the threshold enhancement effect in decays as B — > BB'M is due to the interactions in the 
baryon-antibaryon system. We will return on this hypothesis in Section III. Rescattering effect have been studied, 
in B decays as well as in J/ip decays, in some recent papers based on potential-like model for baryon-antibaryon 
interaction (Paris or Bonn potential) or by using one-pion-exchange model (OPE) [22-24]. In this paper we address 
the final state interactions in the framework of Regge theory, in particular we will study in detail the channel B into 
ppir. 



II. FINAL STATE INTERACTIONS FOR THREE-BODY DECAYS 

In this section we write some basic formula for final state interactions model, in the formalism of Regge theory [25], 
for three-body decays. In particular we study in detail B — > ppn. We start to define the kinematic variables as 



{Pp+Ppf 



= (p P + q) 2 , 



v = (pp + q) 2 , 



(2) 



respectively the invariant mass of proton-antiproton, proton-pion and antiproton-pion system; they satisfy the con- 
straint s + x + y = tub + 2m p + . In the following we neglect the pion mass respect to the other ones and we work 




FIG. 2: Final state interactions for three-body decay. Definition of kinematic variables. 



3 



S = 1 + P 

P 



n n n , n,p 

P P,A 
p + 2 p 1 — 



n, p 

P P 



p, n 



FIG. 3: Diagrams that contribute to the rescattering. V is the Pomeron and n,p,A are Regge trajectories. 

with s and x as independent variables. The relation between the 'full' amplitude A, where the final state interactions 
are take into account, and the 'bare' amplitude A is given by the Watson-Migdal theorem [26]: 

A = VSA , (3) 

where S is the rescattering strong interaction S'-matrix for a given partial wave (J = in the present case). In the 
following we consider rescattering without flavour changing, that is our intermediate state is the same of the outgoing 
one. Moreover, at leading order, for completeness one should consider the interaction between particles in all possible 
way as shown in Fig. 2. The S'-matrix (with the initial state = final state) can be written in the following way [27] 



,/A(s,m 2 m 2 ) J\(x,m 2 ml) J A(y,m 2 ,m 2 ) 

5 = 1 + 2i-± T p P(s) + 2i+ T p7r (x) + 2i+ T p7r (y) , (4) 

s x y 

where A(s, m 2 , m 2 ) is the Kallen function and the T-matrices contains the elastic and also inelastic rescattering effects. 
From symmetry considerations we have T p7r — T pv . The contributions to the T-matrices are shown in Fig. 3. 

The leading contribution to the T-matrices, the Pomeron, can be parameterized, for small and negative transferred 
momentum t, in a universal form as [281 



V(s,t) = -P{t) ij-) e -*W«, (5) 

with so « 0(1 GeV 2 ) and a-p(t) = 1.08 + 0.25 1; this form is suggested by fit to hadron-hadron scattering total 
cross sections. The prefactor (3{t) in Equation (5) represents the Pomeron residue, that is Pomeron 'coupling' to the 
hadrons; phenomenologically it is given by [25, 29] 

The coefficient f3 v writes in term of 0^ u , the V— quark coupling (we consider SU(2) symmetry), and applying the 
additive quark counting rule [28] at each vertex; for instance for the first and the second diagram in Fig. 3 we have 
respectively: 

Ppp = (3 C) 2 . P% = (3 C) (2 C) - (7) 

where 0% u = 1.87 GeV" 1 [29]. 

The non-leading contributions to the rescattering are given by Regge trajectories. These amplitudes can be param- 
eterized by means the Veneziano formula [30, 31] in the following way 

1 -4- (— 1 

TZ(s, t) = -p n +l > 2 6 r (In - a n (t)) (a') 1 ^ (a's)^ , (8) 

where T is the Euler function and 

an(t) = s n + a'(t - m\) , (9) 

is the Regge trajectory with the universal slope a' = 0.93 GeV -2 . In Table I we list the parameters s-n, and m-jz 
for the trajectories n, p and A. f3 n in Eq.(8) is the residue for the Regge trajectory and, as in the case of Pomeron, it 
factorizes in a product of two residue one for each diagram vertex. To estimate the Regge residue we note that near 
t ps m 2 ^, Eq.(8) reduces to a Feynman-like amplitude [33] 

ft(M) « -^Xp , (10) 
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TABLE I: Parameters of the Regge trajectories. 



which allows us to identify f3 n as a product of coupling constants, one for each vertex. 

To estimate the A-proton-pion coupling we used an effective lagrangian approach with the experimental value of 
A — > pn decay width. The interaction lagrangian can be written as 

-^ApTT = gApTT-ip Xfi <9 M 7T + C.C. , (11) 

where -0 is the standard fermionic field for the proton, tt the pion field and x M the Rarita-Schwinger tensor [34] for a 
|-spin field that describes the A. From Lagrangian in Eq.(ll) one can calculate the following decay width 

r(A — > pir) = — P ) g l p7r , , (12) 

and by means his experimental value, T exp (A — * jot) = 118 MeV at p/ = 229 MeV [32], one can determinate the 
coupling constant 

5A^ = 15.4 GcV- 1 . (13) 

In the following we will neglect the i-dependence of (3^ Regge residue and we will identify it as <7A p7r . The other 
residues useful in our treatment are given by [31] 

/3; p = 25.2, /3£ p = 13.0, ^ = -8.0. (14) 

Therefore the T-matrices present in Eq.(4) write as [29, 35] 

T-» = ±- 3^^/_° s+4ro? ^(M) + dt (15) 

^(x) - 1— / (m2 _ x)2 ^(x,i)+^(x,i)+< T (x,i)^, (16) 

lb7r A(x, nip, m^.) y ( p ' 

and finally the double differential width for B — > 7rpp, with the final state interaction effect take into account, is thus 
given by [35] 

1 1 ■ .,2 



dT ' ^ (2^32^^ ^ 

\S\\A\ 2 dxds. (17) 



1 1 ~ M2 



(2tt) 3 32m| 1 



III. DISCUSSION AND RESULTS 



The double differential width writes in Eq.(17) depends on two functions: (i) The bare amplitude A that describes 
the B decay vertex and (ii) the S*-matrix S that describes the re-interaction between the outgoing particles. In 
Section II we discussed the rescattering contributions; the problem now is to understand the role that plays the bare 
amplitude in the full dynamics. How we saw in the Introduction experimental data show that the characteristic peak 
at the threshold in the proton-antiproton system as well as in other different baryonic systems has a very universal 
beaviour. On the other hand the bare amplitude should be strongly dependent from the dynamics of weak decay 
vertex, or in other words, it should be different for different final states. While the strong interactions are independent 
from dynamics of particle productions and they are very similar beaviour between hadrons. Therefore, how previously 
argued in the Introduction, one should realize that the dynamics of the process is principally dominated by rescattering 
interactions between baryons and B decay vertex plays a 'fine-tuning' role on the shape of differential width. 
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FIG. 4: Differential width of B — > ppir as function of invariant mass of proton-antiproton. Shared yellow region is just the 
phase-space $(s). Blue line refers to the final state interactions effect with only re-interaction in the proton-antiproton system 
(T pp ), while red line take into account full rescattering (T pp + 2T pn ). Units are arbitrary. 
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FIG. 5: Dalitz plot of B — > ppix with re-interaction effects. We can distinguish two regions in the diagram: First on the left due 
to the Pomeron exchanged in proton-antiproton system and one at bottom on the right that is the exchange of the A between 
pion and (anti)proton. 



In summary we neglect the variations of bare amplitude respect to S-matrix and we trait A as was a constant. Thus 
if we turn-off the rescattering effects (S — 1) and we integrate-out the x variable in Eq.(17), one obtains a differential 
width proportional to three-body phase-space namely 

$(«) = 5i 3 3 m VV a ~ 4m P' (18) 

that the bare amplitude can only modulate. 

The effects of final state interactions on the differential width of B — > ppir (modulo \A\ 2 ) are reported in Fig. 4 
where the blue line takes into account only the effects of re-interaction between proton-antiproton (see Eq.(15)) while 
red curve contains the rescattering with the pion too (Eqq.(15) and (16)). There are two characteristics that emerge 
from full rescattering in Fig. 4: (i) The strong peak near to the threshold area and (ii) the 'see-saw' trend for higher 
values of s. Both are present in experimental data, compare for instance Fig. 1. These features are better understood 
if we look at the Dalitz plot in Fig. 5; how one can see from the structure of the Dalitz distribution, respect to the s 
variable, the peak is essentially due to the Pomeron for rescattering between the baryonic pair (structure on the left 
of plot) while the latter is due to the presence of the Pomeron and A in the system pion-(anti)proton (structure on 
the bottom). 

To further test our analysis we applied rescattering approach to the channels B — > ppK, B — > ppK* and B — ► ppD. 
In these cases for the systems pK, pK* and pD we take into account only the Pomeron contribution; results are in 
Fig. 6. How one can see also in these cases the rescattering model is able to reproduce the peak near to the threshold 
area but in some cases, as in the channel B — > ppK* for instance, the spectrum is not well reproduced for higher s: 
Probably in this case one should add another rescattering contributions between K* and (anti)proton, but the peak 
is, in any case, due to the rescattering between baryonic pair. Moreover the shape of spectrum is similar in the three 




FIG. 6: Differential width (left to right) of B —> ppK, B — * ppK* and B — + ppD as function of invariant mass of proton- 
antiproton. Shared yellow region is just the phase-space. Rescattering contributions in proton-antiproton system are the same 
of Fig. 4 (blue line), while for full rescattering (red line) the systems pK, pK* and pD contain only Pomeron contribution. 
Units are arbitrary. 




FIG. 7: Differential width (left to right) of B — > J/tppp and B — > J/ifiAp as function of invariant mass of baryon-antibaryon. 
Shared yellow region is just the phase-space. Rescattering contributions are only in baryon-antibaryon system. Units are 
arbitrary. 



cases, but as the phase-space of s variable become smaller (i.e. the outgoing mass of meson increase) the peak width 
tend to enlarge. 

There is another support to our thesis of the re-interaction, and it comes again from experimental data. In fact 
recently BELLE studied B decay channels with an heavy particle in the final state, precisely B — ► J/tpAp and 
B — > J/ippp [36]. From these measurements emerged that there is not threshold enhancement effect at all in the 
invariant mass of Ap when there is a heavy meson in the final state (on the contrary the peak return if meson is 
light one as in the case B — > Apir). Our results for B — > J/ipAp and B — > J/ippp are reported in Fig. 7: For these 
cases we consider only the rescattering between baryon-antibaryon system because the phase-space of x variable is 
small and it gives a negligible contribution to the full amplitude; moreover in the system antiproton-A we consider 
the contribution of Pomeron and the Regge trajectory of p. These two channels are the limit cases of B — > ppD where 
the phase-space of s is so small that only the peak is present. 



IV. CONCLUSIONS 



In this paper we presented a simple model to understand the threshold enhancement effect present in the spectrum 
of many baryonic three-body B (and J/ip) decays. This model is based on hypothesis of re-interactions between 
outgoing particle. We studied in detail the B — > ppn channel because in this case we are able to make a full analysis 
of main rescattering diagrams in the proton-antiproton system as well as in the pion-(anti)proton system, but the 
results can be easily extend to other decay channels. Our principal ansatz is the Constance of bare amplitude, or 
better, its variations from a channel decay to another are negligible respect to the rescattering effects. This feature is 
suggest from experimental behaviour. 

In this framework we are able to reproduce the experimental enhancement effect near to the threshold area for 
B — ► ppir decay (see Fig. 4) and with some change also in other decays where the pion is substituted by a heavier 
meson (see Fig. 6). We studied also two channels with a heavy meson in the final state, precisely B — > J/ippp and 
B — > J/ipAp (see Fig. 7): in these cases due to the phase-space configuration the peak disappear, or better it span all 
spectrum. This feature has been experimentally observed in B — > J/ipAp decay; on the other hand, for B — » J/ippp 
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it is a theoretical prediction. 
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